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Kinetic difTusion of cosmic rays ahead of perpendicular shocks induces charge non-ncutrahty, 
which is mostly, yet not completely, screened by the bulk plasma via polarization drift current. 
Hydrodynamic shear instabilities as well as modified Buneman instability of the polarization current 
generate the turbulence necessary for a Fermi-type acceleration. Thus, similar to the case of parallel 
shocks, in perpendicular shocks the diflfusing cosmic rays generate unstable plasma currents that in 
turn excite turbulence. This allows a self-consistent evolution of a shock-cosmic rays system. In 
the kinetic regime of the modified Buneman instability, electrons may be heated in the cosmic ray 
precursor. 

I. INTRODUCTION 

Acceleration of cosmic rays is one of the main problems of high energy astrophysics. Shock acceleration is the 

leading model [3]. Particle acceleration at quasi-parallel shocks (when the magnetic field in the upstream medium 
is nearly aligned with the shock normal) and quasi-perpendicular shocks (when the magnetic field in the upstream 
medium is nearly orthogonal to the shock normal) proceeds substantially differently. Most astrophysical shocks are 
quasi-perpendicular, yet theoretically acceleration at this type of shocks is less understood than in the case of quasi- 
parallel shocks. It is recognized that the feedback of accelerated cosmic rays may considerably modify the parallel 
shock structure [1]. 

For parallel shocks, the key issues is the generation of magnetic turbulence in the upstream plasma, required to 
overcome the particle escape ahead of the shock. A powerful MHD instability is driven by cosmic rays themselves 

(Bell's instability [2]), providing a self-consistent description of the cosmic ray acceleration, in a sense that that the 
turbulence and accelerated particle may reach a steady state. (The transitional period of reaching that state, i.e. 
starting acceleration from first principles and particle injection, is not addressed by these models). 

In case of perpendicular (super luminal) shocks, the main problem is advection of particles downstream away from 
the shock: in such shocks particles cannot catch up with the shock by streaming along the field lines. Here diffusion 
in the downstream region is assumed to play a major role. Mostly discussed is a particle diffusion due to field line 
wandering [8, 12]. 

Here we investigate a model similar in spirit to Bell's [2] approach, but in applications to perpendicular shocks: 
assuming that a shock accelerates particles that diffuse kinetically {not by field line wandering) far upstream, we inves- 
tigate the consequence of this assumption. The key difference between cosmic ray-modified parallel and perpendicular 
shocks is that the cosmic ray precursor in perpendicular shocks is non-neutral. A shock picks up and accelerate some 
fraction of ISM ions; thus some ions, which were behind the shock if diffusion were absent, are transported ahead of 
the shock. This creates excess of positive charges upstream and a deficit downstream. We stress that we investigate 
a possibility that ions diffuse many Larmor radii upstream; this is different from the problem of electron heating and 
acceleration by reflected ions, which occurs on one Larmor radius of reflected ions [5, 14],. 

In case of parallel shocks such excess of positive charges upstream is balanced by the electrons drawn along the mag- 
netic field into upstream plasma, driving Bell's [2] instability, that excites plasma turbulence and provides scattering 
centers necessary for Fermi-I acceleration. In contrast, in case of perpendicular shocks, the charge disbalance cannot 
be easily cancelled, since kinetic diffusion of electrons across magnetic field is much weaker than that of high energy 
cosmic rays. In principle, the charge density induced by cosmic rays can be canceled by electrons coming "from the 
sides" . But in the frame of the incoming plasma, the charge density builds on very short time scale of the order of 
diffusion length L over shock speed Vg, so that for sufficiently large transverse dimensions of a shock lj_ > VT,eL/vs, 
electrons transported along the magnetic field at electron thermal velocity vt.c does not have sufficient time to cancel 
the ion charge. Thus, even non-relativistic shocks are expected to develop charge density upstream, which, as we 
show below, play an important role in determining the structure of the cosmic ray precursor and possibly leads to 
turbulence generation. 

Cosmic rays diffusing ahead of the shock offset the charge balance in the incoming plasma, which becomes non- 
neutral, with electric field directed along the shock normal. The incoming upstream plasma will partially compensate 
this charge density by a combination of electric and polarization drifts. First, the electric field created by cosmic rays 
will produce electric drift along the shock surface. All plasma components will drift with the same velocity, so that 
electron and ion currents cancel each other, while the cosmic ray current contributes to enhancement of the magnetic 
field in the frame of the shock; magnetic field in the fiuid frame remains the same, see §111 B. Since the electric field 
increases towards the shock, while the rest frame magnetic field remains nearly constant, the electric drift velocity 
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increases. This bulk acceleration drives an inertial polarization current along the shock normal, whereby ions slightly 
accelerate toward the shock, their density decreases, partially compensating the charge density induced by cosmic 
rays. Thus, the charge disbalance created by cosmic rays is partially compensated by polarization current in the 
upstream plasma. 

There is a number of instabilities that can operate in the cosmic ray-modified upstream plasma. First, there are 
fluid instabilities due to shear in the upstream plasma. If the flow is inviscid, the Ray leigh- Taylor instability due to 
variable electric fleld drift may be excited if there is an inflection point in the electric drift velocity (Rayleigh criterium). 
Secondly, small viscosity may drive viscose instabilities, if the Reynolds number is high enough. Finally, current-driven 
instabilities, in particular of the modified Buneman type, may generate plasma turbulence. (We consider the latter 
in more details below.) In all these cases, the resulting instability has wave vector preferentially perpendicular to the 
initial magnetic field, generating the field line wandering required for acceleration of cosmic rays in the first place. 
Thus, similar to parallel shocks, assumption of turbulence and cosmic ray acceleration leads to turbulence generation 
by cosmic rays themselves. This indicates that a steady state of cosmic ray-modified shocks is possible. 

II. PRINCIPAL ISSUES 

A. Distribution of cosmic rays ahead of the shock 

Let us neglect, as a first approximation, the dynamical infiuence of the cosmic rays on plasma bulk motion. In the 
frame of the shock, plasma with density no,/SM is moving with constant velocity Vg in the positive z direction. Let the 
shock be located at 2 = 0, see Fig. 1. Consider cosmic rays in the fiuid approximation, described by a local density 
Tier (this implies that we consider scales much larger than Larmor radii of cosmic rays). Cosmic rays are advected 
in the positive z direction with velocity Vg and diffuse with a given diffusion coefficeint k. Balancing advection and 
diffusion, 

VsdzUcr = Kd^ricr (1) 

we find 

ricr = no,cre''^^, z<0 (2) 

where L = n/vg and we assumed that the shock accelerates cosmic ray protons with typical density no,,.,-- A total 
excess positive charge of cosmic rays per unit area ahead of the shock, a = no,crL is compensated by the similar lack 
of positive charges behind the shock. As a first approximation we neglect changes in the diffusion coefficeint k due 
changing upstream magnetic field, so that L is a given parameter of the problem. As an estimate (an upper limit, in 
fact) we can use Bohm approximation for the diffusion coefficient, k = 'yc'^/u)B,i, where tOB,p = eB/{mpc) is proton 
cyclotron frequency and 7 is a typical Lorentz factor of cosmic rays. In this case the diffusion length 

L ^ jc^/{ujB,iVs) (3) 

is larger than ion Larmor radius in the shocked plasma tl ~ Vs/oJB,i by a factor '){c/vsY ^ 1. Also, diffusion length 
is much larger than cosmic ray Larmor radius by a factor c/vg » 1. 

Previously, a number of authors discussed back reaction of ions reflected from the shock into upstream medium 
[e.g. 13, 17]. Reflected ions propagated back into the upstream plasma a distance of the order of one Larmor radius 
tl ~ Vs/'jJB,i [6]. This may be important for dynamics of electrons {e.g., pre-acceleration) , having small Larmor radii 
[5], but the reflected ions do not affect plasma motion on scales larger than Larmor radius of plasma ions. In contrast, 
we assume that cosmic rays can kinetically diffuse on much large scales n, and investigate the consequences of 
this assumption. 

B. Heuristic derivation of upstream dynamics 

As a result of cosmic ray diffusion ahead of the shock, an electric field is created in the z direction. This electric 
field will induce a drift velocity along the shock plane, which in turn will induce polarization drift that will partially 
cancel the initial electric field. Let the resulting total charge density be ritot and corresponding electric field given 
by AbfEz = Aireritot- Let us neglect variations of magnetic field in front of the shock. (A more detailed derivation 
is given in §111 C.) The resulting electric field and the initial magnetic field induce a drift velocity in the positive y 
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direction Vd = {Ez/ B^jcGy, sec Fig. 1. In addition, plasma components will experience polarization drift, which is 
proportional to mass and thus is mostly important for ions; 

Vp = ^ Vsd;,E = 'inentot = (4) 

Continuity of ion flow requires 

ni{vs + Vp) = nojsMVS, (5) 

thus creating an excess ion charge 

V • 
Sni = — -riQjsM = -ntot-j^ (6) 



At any point the total charge is 



. n-cr ^B,% 
ntot = - OUi = J— « ricr^— < Ucr (7) 

1 H — ¥^ P'* 

Thus, the polarization current cancels most of tlw^ charge induced by cosmic rays. 
Using Eq. (7) we find the velocity of the polarization drift 

Vp Vg^cr 2 I 2 ^ ^sCcr (^) 

'^p,i + ^B,i 

Thus, if a fraction of accelerated cosmic rays ^cr is small, the polarization drift is much smaller than the shock velocity. 
In case of exponential distribution of cosmic rays, the resulting electric field is 

« 47reL'^ner(2) (9) 

and the resulting electric drift is 

Ve ^ ^cr{z)LwB,i (10) 

If the diffusion length is given by Eq. (3), 

Ve ~ icr{z)"f— (11) 
Vs 

We should also verify that the polarization drift is not affecting strongly the cosmic ray distribution. For the density 

(7) the polarization drift of cosmic rays is 

Vp,cr « £,cr{z)jVs (12) 

It can be neglected if ^cr < l/l- 

The above relations give simple estimates of the effects induced by cosmic ray diffusion ahead of the shock. The 
more detailed derivation is given in the next Section. 

III. DRIFT MOTIONS IN NON-NEUTRAL PLASMA 
A. General relations 

Let us consider stationary plasma motion in the frame of the shock. Let ion, electron velocities, magnetic and 
electric fields be given by 

Vj = {0,Wi,j„Wi,J 

Ve = {0,Ve,y,Ve,z} 
E= {0,-Eyfl,E,} 

B = {5„0,0} (13) 



L 



A Bo <J'f^using 




Shock 



FIG. 1: Plasma flows in he frame of the shock. Far upstream, the incoming plasma moves with velocity Vs, magnetic field is 
along X direction and inductive electric field along —y direction. At the shock, cosmic rays are accelerated and diffuse ahead of 
the shock a typical distance L, creating an electric field along the shock normal. Electric drift of cosmic rays in this induced 
electric field and the initial magnetic field produces ^-dependent electric drift ve in y direction. Acceleration of plasma in y 
direction in turn produces polarization drift of ions in z direction. 



{Eyfi is constant because curlE — 0). Assuming massless electrons and massive ions, equations of motion are 



E 



X B = 



(v,V)v, = — (e + ^ X B 



(14) 



Maxwell's equations and continuity require 



curlB = I njVj + Ucr \—c ^ ey - rigVe 
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divE = 47re(ni + ricr — ne) 

riiVi^z = nojSMVg = neVe,z (15) 

Cosmic rays experience just the electric drift along y direction due to combined electric field of cosmic rays and plasma 
particles. 



Eliminating electron's velocities from Eq. (14) 



Ey 

Ve,z = 



(16) 



we get four equations for Vi^y,Vi^z, B^, E^: 

e 



Vi,zVi 



y 



rrii 



Vi.zv'i ^ = — {e. 



rrii 

E'^ = A-Ke ( noisMVs ( — ) + ric 













Bx 




CEyfi 




Ez ■ 


\'"i,z 





D/ 47re / / Vi^y Ez \ ^ Ez \ 

B. Zero inertia limit 

In this limit we assume that rrii ^ 0. This limit is typically not applicable to astrophysical plasmas, since it implies 
that u}p^i/oJB,i 0, while, in fact, an opposite limit, 0Jp^i/ijjB,i ^ 1 is relevant. Still, the limit of massless protons 
allows an exact treatment and highlights a number of important details. A somewhat different approach to this 
problem in given in Appendix A. 

In the massless limit, Wj 0, we find 

Ez 

-Vi 



^yfi 



Vi,z 

which gives 

B' 



= ^ (18) 



Bx 

E'^ = Airencr (19) 



For a given distribution of cosmic rays, of the form ricr = rio^cr^^^^^ we find 

1 

~ ^l + {L/rB{z)y 
L/rsiz) 



--Vs 
\2 



Bx = BoVi + WMzW 



^l + {LlrB{z)YEz = -^Bo (20) 



rniz) 



where 
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is the so-called magnetic Debye radius associated with cosmic rays. (The statement in [9] that in non-neutral plasmas 
the magnetic field is screened on the magnetic Debye radius rs is incorrect). Thus, if L ~ rB(0), the electromagnetic 
drift can become relativistic. 

Note that magnetic field in the plasma rest frame B' = B^/T, where F = 1/^1 — [Vy + v1)/c^ remains constant, 

B' = Bo^/T^^JyT/cj^ (it would have been more correct to use B' instead of Bq as a parameter of the problem). 
For a typical ISM plasma 

rs « lO^cm^B.snoJsM (22) 

is very small even for small cosmic ray fraction ^ 1- This highlights the fact that even a small charge disbalance 
can lead to substantial changes in the upstream region. 



C. Drift equation with polarization current 

In the previous section we showed that charge disbalance ahead of the shock can lead to substantial velocities, and 
thus requires taking inertial contributions into account. In Eq. (17), eliminating y component of ion velocity 



Bx eB:r. 



(23) 



we find 



Vi,z -Vs = {b^dzEz - dzbxEz) ' ^ - {{dzVi^zf + Wi.z^^Wi.z - Vt,zdzVi,zdzbx) 2n2h2 

dzbx = ( "()(— - bx) + Tier] - minoVsdzVi, 

dzE, = Aire (no{^ - bx) + n^r^ (24) 

where we introduced bx = Bx/Bq. 

The drift approximation corresponds to expansion in 1/e, one over electric charge [11]. In a non- neutral plasma it 
is important to keep two leading expansion orders: 

Vi,z -Vs = {bxdzEz - dzbxEz) * ''"^ 



eBlbl 



47r f^f , Vs 

Q^x \ \ ^i,z 



dzbx = -^2^ ( Ez [ noi- bx) + Ucr ] - minoVgdzVi, 



V,, 



dzEz = 47re ( no(— - bx) + ricr ) (25) 



To proceed further, we assume that the polarization drift is weak, and expand in small quantities Vi^z — Vg, bx — 1, 
Ucr and Ez- 



Vi,z -Vs = -Vs{bx - 1) + -^^dzE, 



e 



-'z-^z 



Which immediately gives 



dzbx = ^ OzVi^z 

dzEz = Ane I no{— — —-bx) + ncr] (26) 



^^^^^^ = -Mf{bx - 1) (27) 
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where Ma = Vs/va- The solutions for the electromagnetic fields and components' velocities are 

■ 1 , icr{z) 



(l + /3i)(l-l/Ml) 

Lu)B,i ^criz) LU)B-, 



J- + Pa 



Lu>B,i 1 + 2 T1 ^cr{z) 



ne = riojSM \^ + 
Tli = noTSM I 1 + 



(l + /3^)(l-l/Ml) 

(l-l/Ml)(l + /3i) 
(Ml + l)(l + /3i) 

'^tot = 1 I ^^2 "cr ~ /^A'^cr < ricr (28) 

i + Pa 

where /3a = iijJB,i/i^p,i)- (Relations (28) assume a specific cosmic ray density distribution, ^criz) = ^cr,oe^^^.) 
Equations (28) solve a problem of modification of the upstream plasma in perpendicular shocks by cosmic rays 
diffusing ahead of the shock. 

The electric drift velocity is typically much smaller than the shock velocity. If we introduce parameter r] as cosmic 
ray acceleration efficiency, 

'yncr'tniC? 

r]= 2/0 ' 29 

then with partial screening of the cosmic ray-induced electric field by the incoming plasma, the transverse velocity of 
the incoming plasma, 

Ve'^UtC^/Vs'^Vs^Vs (30) 

is typically smaller than the shock velocity. 



IV. INSTABILITIES 



A. Fluid instabilities in upstream plasma 

There is a number of instabilities that can operate in the upstream plasma. First, there are fluid instabilities 
that may develop in the sheared flow in the upstream medium: KelvinHelmholtz (KH) instability of ideal flows and 
viscously-driven instabilities of flows with high Reynolds numbers. KH instability of sheared flows requires inflection 
point in the velocity profile (Rayleigh condition). This cannot be achieved in the present model, where diffusion 
coefficient was assumed to constant, so that electric drift velocity is proportional to smoothly decreasing density of 
cosmic rays, Eq. (28). On the other hand, if diffusion has a more complicated spacial dependence {e.g. it is higher 
closer to the shock, where the level of turbulence is higher, one might expect an inflection point. 

Allowing for spacial variation of the diffusion coefficient, the cosmic ray number density becomes 

din = ^^^^5,n (31) 

Thus, if the diffusion changes on scale L ~ v/k, one may have an inflection point in the flow. (Note that since diffusion 
is lower far away from a shock, dzK > 0.) For example, if the diffusion coefficient has a dependence k = kqc^^^, the 
inflection point occurs at z = —Lln{Ko/{Lvs)). Eq. (28) then indicates that the electric drift velocity will also have 
an inflection point, resulting in KH instability. 
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B. Electromagnetic instabilities: modified Buneman 

In addition to fluid instabilities, plasmas with "real" drift of particles (as opposed, e.g. to effective Larmor drift 

in inhomogcncous plasma), is subject to a number of powerful instabilities, related, generally speaking to Buneman 
current instability, and in particular to the so called modified Buneman instability [4]. 

In keeping with our approach, we next consider hydrodynamic-type (non-kinetic) instability driven by a perpen- 
dicular ion current. Since the modified Buneman instability is of the mixed electrostatic-electromagnetic type, it is 
necessary to keep electron inertia. Let us transform to a frame which is drifting with velocity ve along the shock 
normal and is advected with shock velocity Vs- In this frame the ions are drifting with respect to electrons with the 
polarization drift velocity Vp ~ ^cr^^s directed along z axis. In considering stability of such plasma, we neglect small 
variations of the magnetic field and associated current along y direction, as well as related variations of electron and 
ion densities. 

Since the electrostatic contributions to the dielectric tensor are most important, we chose wave vector of pertur- 
bations along z direction, k = kzBz- Let the fluctuating part of the vector potential (in radiation gauge) A be 

A = ayGy + a^e^. The fluctuating electric and magnetic flelds are (5E = ?'((j;/c){0, fly. a^}. t^B = — ifc^jaj,, 0, 0}. 
Introducing electron and ion displacements = {^e,xiie,vi^e,z} and = {ii,x)£,i,y)ii,z}i the equations of motions for 
electrons and ions are 



9tVe -I- (Ve • V)Ve 

dt^i + (vi • V)vi 

Vi = VpGz + dtii 
Ve = dtie 



(E + Ve X B) 

(E + ViX B) 



(32) 



Electromagnetic flelds are related by 



dtB = -V X E 

dt'E = V xB {mvi - rieVe) 

c 



(33) 



Finally, continuity equations give 



dtUe + V(neVe) = 

dtUi + V(niVi) = 



(34) 



We flnd 



.el,. , 

?e,j/ = 2 2^ y'^^°'y + ^B,eaz) 



^e,z = — 2~ {^B,eay - icoa^) 

nieC — g 

e 1 

■'ba 



crriiijj ijp- — oj- 



{u!B,iay + iuiaz) 



B,i 



u) = u) — kzVp 
Ue = no (1 - ikz£,e,z) 
, , .kzOJ^i,z 

ni = no{l-i 



CO 



(35) 



which gives the equation for normal modes 

n^-l + Se+ ^Si 
-iD 



iD 

-1 + Se + Si J \a 



(36) 
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where 



Si 



LO^ — UJ7 



D = -^S,+ ^Si (37) 



The dispersion relation follows from (36): 



,~,2 £,2 



The eigenmodes satisfy 



n^ = l-5.-^5.-^^-^ (38) 



— Se — Si , , 

I (39) 



Thus, the eigenmodes are elliptically polarized. Note, that the wave vector of the perturbations lies in the polarization 
plane, so that the wave is of a mixed electrostatic-electromagnetic type. 

For vanishing drift velocity, Vp = 0, in the low frequency limit u)B,e — * oo, Eq. (36) gives the extraordinary mode 



n 



'-^-Se-s.- — ^ = 1 + ^ ^ - ^ y -i,^ . (40) 



In case of finite drift velocity, keeping terms linear in rrie/mi and assuming Co <C ws^j, we find 



n2 = l+(^ + 5^')4^ (41) 



m,; Lip- I UJ 



B,i 



which gives 

a; « fc, (vp ± ^vl-vlirrie/mi + Pl)'^ (42) 
Equation (41) has complex roots provided 



va 



(43) 



If (3 A <C \Jnielnii (this requires S < 7 x 10 ^Gn^/gj^^), the instability occurs for 

Vp > ^/mJrrieVA (44) 
Since polarization drifts velocity is Vp ^ ^cr^s, it is required for instability that 

U > \[^^ (45) 

where AIa = Vs/va S> 1 is the shock Mach number. 

For velocities larger than (44), the instability growth rate is 

/ 777 

r=J— Ccfc.v, (46) 

V mi 

The instability has time to grow ~ L/Vs', the requirement TL/Vg > 1 limits the wave number of growing modes to 

2^mi/me^^ < kJi < 1 (47) 
Mat] 
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where 5 = c/utp^i is ion skin depth. The upper limit on kz < l/6i follows from neglect of kinetic effects. 
Thus, if cosmic ray acceleration efficiency is 

V>'2^f^^ (48) 

instability will have enough time to grow. Condition (48) favors strong shocks with weak magnetic field. 

What is the nature of the unstable modes? Eq. (42) indicates that both eigenmodes become unstable, implying that 
both components of the electric fields orthogonal to the initial magnetic field, as well as fluctuations of magnetic field 
along the initial magnetic field, grow exponentially. Eq. (39) indicates that the phases of the electric field components 
Ez and Ey are shifted by ninety degrees, implying that growing modes correspond to vortical perturbations. These 
kinds of perturbations will shuffle field line along z direction. If we allow for finite perturbations along x and y 
direction, the instability will induce field line wandering. Since field line wandering can bring both signs of charge, 
this might be break down our assumption of charge non-neutral plasma. Still, the motion of electrons and ions along 
the field line might proceed in different regime: resulting magnetic field amplification and corresponding magnetic 
bottles will reflect electrons more efficiently due to their small Larmor radius. 



C. Applicability and saturation 

The applicability of our approach requires that (i) the diffusion length L is larger than the Larmor radius of cosmic 
rays rL,cr = 7c/u's,i - this condition is satisfled for any non-relativistic shock, since L ^ {c/vs)rL,cr ^ i~L,cr\ (ii) L is 
larger than skin-depth of cosmic ray component, 6cr = c/up^cr- We flnd 



Scr VsVA V 2 VA 

This puts a lower limit on cosmic ray efficiency 

r? >2(/?^)2 = 10-10 — B_6' (50) 
nisM 

which can be easily satisfled. 

The growth rate (46) was calculated assuming cold ions and electrons. This is justifled if the polarization drift 
velocity is larger than electrons thermal speed vt^bj Vp = ^cr^s > ■\/T/me. This requires 

(51) 

where Mg = Vs/vT,i is the conventional Mach number. The condition (51) is quite restrictive, so that the modifled 
Buneman instability may also heat electrons. Previously, it was suggested [7, 15] that lower hybrid waves in the 
CR precursor of a perpendicular shock might be a plausible electron heating mechanism. To assess the efficiency 
of electron heating, kinetic-type calculations are required; we postpone them to a future paper. For drift velocities 
smaller than ion thermal speed, Vp < VT.i it is expected that the instability will be stabilized. This requirement places 
a lower limit on cosmic ray faction, ^cr > ^/Ma- 

To estimate the saturation level of instability, we assume that a large fraction of the ion polarization drift is 
converted into fluctuating magnetic fleld. This gives an estimate 




V. DISCUSSION 



We considered how cosmic ray accelerated at perpendicular shocks and diffusing kinetically ahead of the shock 
modify the upstream flow. First, the electric charge density induced by cosmic ray diffusion is mostly, but not 
completely, screened by the polarization drift in the upstream plasma. The remaining charge density induces shear 
along the shock plane, so that the shock become oblique (in a sense that the flow velocity is not aligned with the 
shock normal; magnetic fleld remains orthogonal to the shock normal). 
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For fast shocks with high efficiency of cosmic ray acceleration (so that the condition (45) is satisfied), modified 
Buneman instability develops. In addition, a sheared flow can be unstable to fluid shearing instabilities. Both fluid 
and drift current instabilities would generate turbulence and presumably enhanced particle diffusion, vindicating the 
assumption that perpendicular shocks are efficient cosmic ray accelerators. The proposed instability can also be 
important for electron heating in the cosmic ray precursor. This requires a kinetic treatment of the instability. In 
addition, since the higher energy particles are diffusing faster, the resulting population inversion ahead of the shock 
may drive cyclotron instabilities; we leave consideration of this possibility to a future paper. 

In case of parallel shocks, development of Bell's instability leads to generation of fluctuating perpendicular magnetic 
fields, which become dominant over the initial parallel field [16]. In that case, the instability discussed in this paper 
can be regarded as a secondary instability of the cosmic ray-modified parallel shocks. 

The effect discussed in this paper is due to non-neutrality of the upstream plasma and thus is easily missed in 
a conventional MHD treatment. In addition, we assume that kinetic cosmic ray diffusion dominates over field line 
wandering, at least upstream of the shock. If particle were to diffuse mostly by field line wandering, charge neutrality 
would be established by drawing a parallel current. Also, we did not address the issues of particle injections, i.e., how 
to start the process going. 

I would like to thank Elena Amato, Benjamin Chandran, Martin Laming and Anatoly Spitkovsky for the most 
enlightening discussions. 
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APPENDIX A: MAGNETIC FIELD IN THE UPSTREAM PLASMA 

In this Appendix we re-derive behavior of magnetic field ahead of the shock neglecting inertial effects. Let us 
assume that a flow carrying magnetic field Bq approaches a perpendicular shock, and that due to the acceleration of 
cosmic ray particles there is a layer of uncompensated charge in front of the shock. Neglecting diamagnetic currents 
induced by cosmic rays as well as pressure of cosmic rays, the equation of motion of cosmic rays and electromagnetic 



12 



fields reduces to 



E + - X B = (Al) 

c 



here E and v are total electric fields and drift velocities. Let the density within the charged layer be ritot and the 
thickness of the layer be L. Since the charge-separated current is j = eritot^-, and V x B = (47r/c)j, we find 



E+ ^- VxBxB = (A2) 



Taking divergence of this equation gives 



c? /V X B \ 

ntot + T^^V • X B = (A3) 

167r^e^ V ntot J 

The Eq. (A3) reminds of the plasma dynamics in the limit of electron MHD [10] . For given ntot and L, the modification 
of magnetic field depends only on number density of cosmic rays, Eq. (A3), and not on their energy. 
Eq. (A3) can be integrated, 

M{z) = [ ntot{z)dz (A4) 

J z 

where C'l and Co arc integration constants. The quantity M{z) is a surface density of cosmic rays ahead of the shock, 
integral of density from a given position z up to the shock front. (Recall that in the zero inertia limit there is no 
compensating polarization drift of bulk ions.) In the particular case when the density of cosmic rays is given by Eq. 
(2), M{z) = no,cr-^(l — e^/^), so that ^/{z = — oo) = no^crL. Since the magnetic field a,t z = — oo is Bq, this gives 

r, 167r^ e^nn rrL „ „ , . , 

Bo = + Cmo^crL + Co A5 

Assuming assume that ^-component of the electric field is zero and magnetic field is Bq a,t z = — oo, we find 

Bl = Bl(l+ 



rsiz) 

= — ^Bo (A6) 

Consistent with Eq. (20). Again, the apparent increase of magnetic field is, in fact, due to the Lorentz boost along y 
direction with electric drift velocity 

VE^ ^ L/rB{z) 
c yi + LVrs(^)2 



